Gold in graphene: in-plane adsorption and diffusion 
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We study the bonding and diffusion of Au in graphene vacancies using density-functional theory. 
Energetics show that Au adsorbs preferably to double vacancies, steadily in-plane with graphene. 
All diffusion barriers for the complex of Au in double vacancy are above 4 eV, whereas the barriers 
for larger vacancies are below 2 eV. Our results support the main results of a recent experiment 
[Gan et a/.. Small 4, 587 (2008)], but suggest that the observed diffusion mechanism is not thermally 
activated, but radiation-enhanced. 
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As the exceptional properties of graphite have been 
long known, it is no surprise that the recent developments 
in graphene fabrication techniques have given an enor- 
mous impetus for graphene research. Not only is pure 
graphene interesting, but also interactions with metal 
atoms, because of their fundamental relevance to appli- 
cations like catalysis, batteries, and nanoelectronics^^i^ 

These interactions have been investigated both exper- 
imentally and theoretically. Previous studies include 
metal atom adsorption and diffusion on surfaces and 
edges of graphite, graphene and carbon nanotubes*^ 
Among experiments, transmission electron microscopy 
(TEM) is a versatile method for sample characterization 
as well as manipulation^iii^, and can achieve nearly single 
atom accuracy. Recent interesting experiment by Gan, 
Sun and Banhart^, observed in-plane adsorption of Au 
and Pt atoms in graphene sheets and carbon nanotubes 
using TEM, and also measured the rate for in-plane metal 
atom diffusion. 

In this Letter we confirm this experimentally observed 
propensity of Au to adsorb in-plane with graphene. The 
experiment suggested thermal in-plane diffusion of Au 
with energy barrier ^ 2.5 eV; our calculations sug- 
gest that the diffusion mechanism is not thermal, but 
radiation-enhanced. 

We used density-functional theory (DFT) with PBE 
functional^^, and projector augmented waves for treat- 
ing the valence electrons^i, as implemented in real-space 
grid code GPAWi^'i^. We started by calculating vacancy 
formation energies, defined for n-atom vacancy as 

£^(form) = [£^(n- vacancy) + n • Scoh] — £^(graphene), 

for single to quadruple vacancies; here Ecoh = 7.9 eV is 
the graphene cohesion energy^i^ Definition above is the 
cost to create an n-atom vacancy and consider the re- 
moved carbons as a part of graphene elsewhere. This 
definition is relevant because it appeared that carbon 
atoms on graphene diffuse rapidly to anneal other vacan- 
cies^^ The vacancy formation energies, shown in Fig. [1] 
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FIG. 1: Formation energies for carbon vacancies in graphene 
(empty squares), and formation energies for Au adsorbed in 
graphene vacancies (filled squares; local geometries shown in 
insets). The difference between the two curves is the Au ad- 
sorption energy. Apart from Au in single vacancy, all struc- 
tures are planar. 



for optimized geometries^^, agree with previous calcula- 
tions For single and double vacancies they are close to 
8 eV, and increase ultimately with ~ 2 eV/C-atom slope 
for larger vacancies. Note that once a single vacancy ex- 
ists, double vacancy costs practically no energy — in nan- 
otubes double vacancy costs actually less than the single 
vacancy/^ The last quadruple vacancy has exceptionally 
low energy due to the absence of dangling bonds in the 
symmetric hole. 

Fig. [1] shows also the formation energies for Au ad- 
sorbed in graphene vancancies, where the reference is 
graphene and Au in vacuum. The difference between 
these two curves, which is the Au adsorption energy, is 
3 — 6 eV; last quadruple vacancy is a clear exception. Sin- 
gle vacancy is too small for Au and optimized geometry 
is not planar, whereas in larger vacancies Au is precisely 
in-plane with graphene. Hence, assuming a planar ad- 
sorption geometry — as demonstrated experimentally — 
and thermal equilibrium, we conclude that Au preferably 



2 



20 




Energy from Fermi level [eV] 

FIG. 2: Projected density of states for s-, p- and d-orbitals of 
Au (above abscissa) and the four neighboring C atoms (below 
abscissa). Insets show the wavefunctions, color standing for 
the phase, for selected s- p- and d-dominated contributions 
for Au. Charging of Au according to Bader analysis is H-0.5e. 



adsorbs in double vacancies in-plane with graphene. 

Why does Au prefer bonding in-plane? This is ana- 
lyzed chemically in Fig. [2l Firstly, the in-plane orbitals 
of Au overlap with the dangling bonds of neighboring 
C-atoms in a bonding configuration; the two p-orbital 
contributions are split due to elongated structure of the 
vacancy. Secondly, and most importantly, out-of-plane 
Au (i-orbitals bind above and below to the 7r-electrons 
of graphene — this gives Au atom stability against out-of- 
plane motion. The in-plane stability was confirmed also 
by molecular dynamics simulations at high temperatures. 
Because the bonding is strong both in-plane and out-of- 
plane, we expect the motion of Au to be confined by large 
energy barriers; in the following we shall investigate Au 
diffusion barriers in a double vacancy in more detail. 

To explore different Au diffusion paths, we used sim- 
ple guessing, we conducted extensive molecular dy- 
namics simulations with density-functional based tight- 
bindingSUiSi, and we used constrained dynamics, for sub- 
sequent analysis. Several paths were optimized, and the 
most prominent ones are shown in Fig. [3l The lowest- 
energy path I with 4.0 eV barrier involves radical out-of- 
plane motion of Au, and path II with 5.8 eV barrier in- 
volves radical out-of-plane motion of C. The lowest nearly 
in-plane path III has much higher barrier of 7.0 eV, path 
IV with barrier even higher. 

We cannot absolutely exclude the existence of other in- 
teresting paths, but we are confident that there are no rel- 
evant paths with ^ 2.5 eV energy barriers, as suggested 
by the experiment.^ Firstly, radical out-of-plane motion 
is unlikely because measured diffusion was quantitatively 
similar for Au in single graphene sheet as well as between 
multiple graphene sheets. Since the barrier is the same 
as Au adsorption energy, overcoming the barrier of path 
I, for example, could equally well result in the escape 
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FIG. 3: Diffusion mechanism paths and their transition state 
energies for Au in double vacancy. Numbers help to visu- 
alize C atom identification during Au jump; gray (blue) for 
C atoms that change positions. Paths were optimized with 
nudged elastic band method ^^'^^ starting from almost linearly 
interpolated initial guess. Path I contains out-of-plane motion 
of Au, whereas path II contains out-of-plane motion of C(#l). 
Path IV contains only translation, other paths also rotation. 



of Au. Secondly, any simple site-exchange mechanism or 
out-of-plane motion of Au is particularly unlikely because 
it would require large local momentum imbalance due to 
large mx^/mc = 16.4 mass ratio. Thirdly, the lowest of 
the nearly in-plane paths has 7.0 eV barrier — more paths 
than shown were examined. Since the in-plane motion of 
Au requires breaking several bonds, and since graphene 
is exceptionally stiff material, nearly three-fold reduction 
in energy barrier, for some unforeseen path, is not plau- 
sible. 

On the other hand, the barrier for a single vacancy 
that merges with double vacancy containing Au is only 
0.2 eV. Since the triple vacancy, as shown in Fig. [H has 
extra space around Au, in-plane motion is easier. And 
indeed, it was not difficult to guess a strictly in-plane 
diffusion path for Au with barrier below 2 eV. Hence, in 
triple vacancies Au could diffuse easily; but Au existing 
in triple vacancy cannot be justified on thermodynamic 
grounds, and, further, the barrier is too low and diffusion 
would have been too fast for TEM observation in the 
experiment of Ref. 0. 

Now we shall juxtapose our calculations with the ex- 
periment by analysing information given by Can et al^ 
To begin with, we point that the barrier 2.5 eV was 
not deduced from Arrhenius-type behavior, but from the 
equation 

D = ga^VQ exp {-Eb/ksT), 
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where D is the diffusion constant, measured by directly 
monitoring atoms, ^ ~ 1 is a geometrical factor, a is the 
lattice constant and vq is the attempt frequency. The 
measured diffusion constant for Au at T = 600 °C was 
D = 6 X 10-22 - 2 X 10-21 m^s-^ Taking a = 1 - 
2 A(depends on mechanism) and ^ 5x 10^^ g-^,^^ the 
time between Au jumps for two-dimensional diffusion is 
- 2 - 20 s. 

How about radiation? With the carbon atom displace- 
ment threshold of 15 eV, the radiation from 30 A/cm^ 
current density was estimated to cause displacement of 
every carbon atom in 180 s, where most of the vacancies 
are rapidly annealed.^'22 Hence, from selected N atoms 
at least one is displaced every ISO/N s. Assuming Au 
in double vacancy to have 14 neighboring C atoms, at 
least one of these is displaced every ~ 10 s. Vacancy's 
environment is energetically excited at this rate, and the 
excitations may help Au to overcome large energy bar- 
riers. This rate, obtained from radiation effects, is too 
close to independently obtained Au diffusion jump rate 
to be a mere coincidence. Our conclusion, based on our 
calculated energetics, barrier heights, and the experimen- 
tal analysis, is that the diffusion in the experiment is 
radiation- enhanced diffusion^ 

The radiation-enhanced interpretation would explain 
some peculiar features in the experiment. First, the dif- 
fusion constant for Au and Pt was similar, which was not 
expected on the account that C-Pt interaction is stronger 
than C-Au^ This can be now understood: if the mech- 
anism is dominated by radiation enhancement, the C- 
metal interaction will not play a crucial role. Second, 
one-dimensional diffusion, and diffusion in the curved 



layers of carbon nanotubes, was very similar to pure 
two-dimensional diffusion. This cannot be understood if 
the diffusion is thermal — it is not plausible that different 
chemical interactions give the same barriers for Au and 
Pt in different environments. Using radiation-enhanced 
argument, in turn, carbon environment around the metal 
can be thought to be excited occasionally into extremely 
vigorous motion, causing insensitivity to interactions and 
geometries. We shall continue our calculations with Pt 
and other transition metals. 

What could the actual diffusion mechanism be? The 
mechanism could be diffusion of Au in double vacancy. 
Radiation can frequently provide the required 4 eV into 
the neighborhood of Au — even though every push cer- 
tainly does not cause a jump. The mechanism could be 
merging of a single vacancy with Au in double vacancy, 
followed by diffusion of Au in triple vacancy, and finally 
followed by annealing of triple vacancy with diffusing C 
atoms. Complex mechanism like this would explain the 
temperature-dependence that was observed in Fig. 4 of 
Ref. M' There are several possibilities for the mechanism, 
once it is radiation-empowered. 

In conclusion, because the problem has technological 
relevance, and as the ultimate solution is hard to provide 
by further calculations, the status calls for experiments 
to clarify the relations between thermal and radiation- 
enhanced diffusion. 
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